typeset using JPSJ.sty <ver.l.Ob> 



Optical Conductivity of the Trellis-Lattice t-J Model: 
Charge Fluctuations in NaV 2 Os 

Satoshi Nishimoto 1 and Yukinori Ohta 1 - 2 

1 Graduate School of Science and Technology, Chiba University, Inage-ku, Chiba 263-8522 
2 Department of Physics, Chiba University, Inage-ku, Chiba 263-8522 

(Received August 6, 1998) 

Optical conductivity of the trellis lattice t-J model at quarter filling is calculated by an exact- 
diagonalization technique on small clusters, whereby the valence state of V ions of NaV^Os 
is considered. We show that the experimental features at ~1 eV, including peak positions, 
presence of shoulders, and anisotropic spectral weight, can be reproduced in reasonable range of 
parameter values, only by assuming that the system is in the charge disproportionated ground 
state. Possible reconciliation with experimental data suggesting the presence of uniform ladders 
at T > T c is discussed. 
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The valence state of V ions in NaV205, a compound 
once believed to be a typical spin-Peierls system,ErEf has 
become an issue of intense controversy even above the 
transition temperature T c ; mpeirLX-ray structural anal- 
ysis and other experimentgj'Q'Eltf seem to suggest the 
uniform valence of V 4 5+ , the system being the coupled 
uniform ladders at quarter filling, whereas an interpreta- 
tion of low-energy optical conductivity data suggests thx? 
system to be in a broken-parity electronic ground state. LP 
To clarify this issue will help one understand the mecha- 
nism of the phase transition of this compound observed 
at T c = 34 K. 

In this letter, we consider this issue by calculating the 
optical conductivity cr(u>) of the trellis-lattice t-J model 
at quarter filling directly. We will show that the experi- 
mental features of <j(oj) at li ~ 0.6 — 2.5 eV, such as the 
positions of the main peak and its shoulders, anisotropy 
in the spectral weight, etc., observedLlEP over a wide tem- 
perature range above and below T c , can be explained 
within a reasonable range of parameter values, only by 
assuming that the system is in the charge disproportion- 
ated ground state. We will thereby point out that the 
difference in the time scale of the measurements may be 
a possible way of reconciliation: charge (or valence) fluc- 
tuation of V ions is slow compared with the frequency 
of electric field of light used in the optical measurement, 
but it is fast enough compared with the time scale of the 
other measurements. We would suggest that measure- 
ments of frequency range of ~50 GHz to ~1 THz should 
be able to detect a possible resonant charge fluctuation 
in this compound. 

The low-energy electronic state of the two-dimensional 
plane of NaV^Os may be written by the trellis-lattice t-J 
model, i.e., coupled anisotropic t-J ladders (see Fig. 1), 
as we have shown in refill The Hamiltonian is given by 

H =- ^(4<^+h.c.) 



Fig. 1. Schematic representation of the trellis lattice t-J model. 
The structure consists of the anisotropic ladders (solid lines) 
connected with zigzag chain bonds (dotted lines). The crystal- 
lographic a- and fo-axes are also shown. 

(ij) 
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where (ij) represents the nearest-neighbor bonds along 
the legs and rungs with the hopping and exchange pa- 
rameters tij and Jij taking t±_ and Jj_ for the rungs and 
f II and J|| for the legs. It also represents the zigzag chain 
bonds with parameters t xy and J xy where the intersite 
Coulomb repulsion V xy is also taken into account. We 
introduce the 'symmetry-breaking' on-site energy of an 
electron, Aj, which takes either A (> 0) or depend- 
ing on sites, whereby we can make the system charge 
disproportionated if A ^ 0. c\ a — c\ a {l — rii- a ) is the 
constrained electron-creation operator at site i and spin 
a ( = tj I); Si is the spin-i operator, and rij = mf + rijj 
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is the electron-number operator. We consider the case 
at quarter filling. 

We employ a Lanczos exact-diagonalization technique 
on small clusters and calculate the frequency-dependent 
conductivity cr(w) to examine the anisotropic optical re- 
sponse of the system. A finite-size cluster of the trellis 
lattice, i.e., coupled two 4x2 ladders shown in Fig. 1, 
is used with periodic boundary condition. The optical 
conductivity is defined by 

c Q (w) = — |(Vv|ja|V'o}| 2 8(w — (E v — E a )), (0.2) 

where \ij) v ) (E v ) denotes the v-\h eigenstate (eigenen- 
ergy) of the system (in particular, v = denotes the 
ground state). Also, j a with a = a, b denotes a compo- 
nent of the current operator 
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ia — i ^ ' tjj (vj rj) a c icr Cj a) 



(0.3) 



where (r, — rj) a is the a-component of the vector con- 
necting between sites i. and j. The system we treat is 
in the insulating stateB' and there is no Drude peak at 
u = 0. 

We use the following two sets of values of the lad- 
der parameters (in units of eVl obtained from the d- 
p model by the method of ref.B> Set A: t± = 0.298, 
t\\ = 0.140, Ji = 0.0487, and J\\ = 0.0293, which are 
obtained in refxP using realistic values of the d-p model 
parameters t p d = 1.22 (1.03) for the hopping integral 
of the rung (leg) V-0 bonds, Ud — 4.0 for the on-site 
Coulomb repulsion on V ions, and A p d — 4.0 (6.5) for the 
energy-level difference between the V d^-orbital and O 
p-orbital on the rung (leg). Set B: t± = 0.397, t\\ = 0.194, 
Jjl = 0.0658, and J|| = 0.0546, which are obtained us- 
ing a different set of values of the d-p model parameters 
with enhanced hopping amplitude: t p d — 1-46 (1-24), 
U d = 5.0, and A pd = 4.0 (6.5) for the rung (leg) V-0 
bonds. Values of the other parameters t xy , V xy , and A 
are then varied for simulating various situations. We use 
the relation J xy — ^t 2 xy /Ud for simplicity. 

The optical conductivity observed in experimen 
has the following features: (i) Onset of charge-transfer 
excitations from O p x - or p y - to V d^-orbitals appears 
at u> ~ 3 — 3.5 eV which is anisotropic between a \\ E 
and b || E. (ii) There appear broad spectral features 
around lu ~ LeV, which have been identified as the d-d 
transitions Jfr LP the main peak has a shoulder in ~0.5 eV 
(~0.8 eV) higher-energy side in a a (uj) (at,(uj)). Also the 
height of the peaks (or the integrated spectral weight) 
is markedly larger in a || E than in 6 || E. (iii) 
There is a small continuum in the low-energy region of 
uj ~ 0.01 — 0.1 eV, the origin of which was argued to 
be 'charged magnons'If The spectral features (i) and 
(ii) are hardly changed over a wide temperature range 
above and below T c Jl3 while the oscillator strength of 
the feature (iii) decreases suddenly at T c with lowering 
temperatures. We now want to discuss the origins of the 
features (ii). 

Let us first take some simplified pictures. Suppose an 
isolated rung of the ladder, then we have the d-d exci- 



Fig. 2. Optical conductivity cr a (w) calculated for the trellis- 
lattice t-J model with A = 0. Solid and dashed curves show 
cr a (u) for a || E and cri,(u>) for b \\ E, respectively. In (a) and (c), 
we use values of the ladder parameters of Set A, and in (b), we 
use those of Set B. All the 5-functions in Eq. (2) are Lorentzian 
broadened with the width of 0.05 eV. Peaks at u) = 0, which 
appear due to the broadening and factor 1/ui in Eq. (2), are 
spurious. 



tation peak in a a (ui) at lu = -J (2t±) 2 + A 2 , the value 
of which is however too small compared to ui ~ 1 eV 
if A = or there is no charge disproportionatiomEP 
because the effective hopping amplitude t± ~ 0.3 eV 
strongly reduced from the bare d-p hybridization t p d is 
relevant here.B' Suppose next the trellis lattice in the 
limit of small hopping amplitude and large V xy and/or 
A where the sites of, say, the right legs of the ladders are 
all with A > and those of the left with A = (result- 
ing in the charge disproportionated ground state), then 
we have a peak in a a {u>) at u> = A + 2V xy - This is be- 
cause the electron transferred through the rung interacts 
(by V X y) with two electrons on the leg of the neighboring 
ladder across the zigzag chain bonds; the peak position 
can be comparable to ui ~ 1 eV because V xy ~ 0.5 eV 
is enough even at A = 0. If small randomness occurs 
in the disproportionated charge configuration, then an- 
other peak can also appear at ui — A + V xy , because the 
transferred electron may interact with only one electron 
on the leg of the neighboring ladder across the zigzag 
chain bonds. Thus, there appears a two-peak structure 
in (7 a (uj) (which, as we will see below, corresponds to 



Optical Conductivity in NaV2 0s 



3 



the experimental broad features around lu ~ 1 eV). It is 
also noted that, if the charge configuration completely 
disproportionated on one side of the legs is assumed, the 
spectral weight for b || E becomes quite small, because 
the electrons are hard to oscillate along the b-direction 
due to the strong on-site Coulomb repulsion (or prohib- 
ited double occupancy in our model) and small values of 
t xy ; any randomness in this charge configuration would 
however make the spectral weight o~b(u) larger. 

Now let us see the calculated results for cr a (oj) with 
keeping these ideas in mind. We first of all find in 
Fig. 2(a), where we use the parameter Set A with A = 
and V xy = (i.e., without charge disproportionation) , 
that the main features appeared at w ~ 0.5 eV are again 
too low in energy compared with experiment. The pa- 
rameter Set B with A = and V xy = also gives similar 
results. We also find that the relative main-peak po- 
sitions between a \\ E and fo || J57 are sensitive to the 
changes in small parameter values t xy and V xy and that 
some broadened features in er a (a->) are induced by the 
presence of t xy . We however find that the essential spec- 
tral features, i.e., the overall main-peak positions, remain 
unchanged. These results thus demonstrate that the 
u) ~ 1 eV features cannot be reproduced in reasonable 
parameter range, without charge disproportionation. 

The charge disproportionated ground state can be 
achieved by including either the large inter-site Coulomb 
repulsion V xy along the zigzag chain bondscP or the on- 
site energy-level difference A between, say, the left and 
right legs of the ladders. In the former case, the cal- 
culated charge correlation functions {niUj) indicate that 
at V X y <; 0.6 — 0.7 eV the strong tendency toward charge 
disproportionation comes outff although the symmetry 
is not really broken in finite-size systems. Figures 2(b) 
and (c) show the results for such situations. Here we 
assume V xy = 1.0 eV, and use values of the parameter 
Set B with t xy = 0.1 eV and Set A with t xy = 0.05 
eV in Figs. 2(b) and (c), respectively. We find in both 
figures that the broad spectral features in CT a ioj) around 
oj ~ 1 eV really appear and some smaller features also 
appear at u> ~ 2 eV. We analyze the final states corre- 
sponding to the largest peaks and find that the features 
at ^1 eV come mainly from the states with 1 lu ~ V xy 
and the features at ~2 eV come mainly from the states 
with 1 lo ~ 2V xy \ as we have discussed above using sim- 
plified pictures. The anisotropy in a a (uj) between o || E 
and b \\ E is also noted; the main peak in (Tb(u) is lo- 
cated in ^0.5 eV higher energy side than that in <J a (ui) 
and the spectral weight for o~b{uj) is considerably smaller 
than that for a a (u>). These results are all in good agree- 
ment with experiment. To get more precise agreement 
with experimental features would necessitates the cal- 
culations in larger size clusters, which are however not 
feasible in cluster diagonalizations. 

Next let us see the results obtained by introducing 
A ^ to break the symmetry and actually realizing 
the charge-disproportionated ground state. Note that 
although such an ad hoc assumption of the presence of 
nonzero A induces charge disproportionation in the sys- 
tem, its origin or the mechanism of appearing A =/= 
should be clarified; the simplest explanationEP may be the 



Fig. 3. Same as Fig. 2 but for the case with A ^ 0. Here we 
assume that the sites on the left legs of the ladders are all with 
A > and those on the right are all with A = 0. Values of the 
ladder parameters of Set A are used. 



presence of strong intersite Coulomb repulsions, which 
induces the broken-symmetry ground state in the ther- 
modynamic limit and leads to A / 0. The calculated 
results for a a (uj) are shown in Fig. 3 where we assume 
that the sites of the right legs of the ladders are all with 
A > and those of the left are all with A = 0. The lad- 
der parameters of Set A are used. We assume A = 0.3, 
V xy = 0.2, and t xy — 0.03 eV in Fig. 3(a), which result in 
the valences of V 4 789+ and V 4 211+ ; here we find that the 
peak positions in a a (uj) are again still too low compared 
with experiment. By increasing the values of A and V xy , 
wc find that the peak positions shift to higher energy side 
and the agreement with experiment becomes better; in 
Fig. 3(b), we assume A = 0.4, V xy — 0.5, and t xy = 0.05 
eV, which result in the valences of V 4 940+ and V 4 060+ , 
and in Fig. 3(c), we assume A = 0.4, V xy = 0.5, and 
t X y = 0.03 eV, which result in the valences of V 4,943+ and 
V 4 ' 057+ . In the latter two cases, we find that (i) the main 
peak appears at u> ~ 1 eV and its shoulders at u> ~ 1.6 
eV in both a a (to) and Ubioj) where the distance between 
the main peaks and shoulders is controlled primarily by 
the value of V xy , (ii) the main-peak position in ct{,(cj) is 
slightly higher in energy than that in a a (ui), and (iii) the 
strong anisotropy of the spectral weight between a \\ E 
and b \\ E appears. These results are all in good quali- 
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tative agreement with experiment. We furthermore find 
that the broad main peak has some internal structures, 
which may correspond to the experimentally observed 
broad peak which also has some structures. t3 

We note that the anisotropy in the spectral weight 
between a \\ E and b \\ E which is seen in Figs. 3(b) 
and (c) is too strong compared with the experimental 
anisotropy. This might suggest that the charge-ordering 
pattern assumed above should be made more 'random', 
so that the electrons may oscillate along the legs of the 
ladders. To check this idea and to improve the strength 
of the anisotropy, we also use ±he charge-ordering pattern 
called the zigzag structurelij'li3 ) which was argued to be 
realistic at T < T c . The calculated result for a a (uj) is 
shown in Fig. 4, where we indeed find that the spectral 
weight of Ub{oj) is much increased without changing other 
features, resulting in a very good agreement with the 
experimental features. 

We have also sought for low-energy (ui < 0.1 eV) fea- 
tures argued as 'charged magnonstP in our calculated re- 
sults. We however could not find any indications: there 
are many low-energy eigenstates but none of them has 
nonzero spectral weight. Interestingly, we find that, if 
we introduce a local lattice displacement, a small weight 
appears at the low energies, which might suggest the rel- 
evance of phonon-assisted processeslr 



other structural analysis experiments, one finds the sys- 
tem to be uniform. It is reported that an ESR exper- 
iment of up to 24 GHz finds no motional narrowing of 
the linewidttO' and that the optical measurement down 
to the lowest frequency of ~10 cm _1 r detects no anoma- 
lous behaviors indicating resonances We would thus 
suggest that measurements of frequency range of ^50 
GHz to ~1 THz should be able to detect the resonant 
charge fluctuation in this material. Precise character- 
ization of the charge fluctuations at T > T c needs to 
be made further from both experimental and theoretical 
sides to clarify the mechanism of the phase transition of 
NaV 2 5 at T = 34 K. 

In summary, we have calculated the optical conductiv- 
ity of the trellis lattice t-J model at quarter filling using 
an exact-diagonalization technique on small clusters and 
considered the valence state of V ions of NaV2 05. We 
have shown that the experimental features observed at 
uj ~ 0.6 — 2.5 eV, including the positions of the main peak 
and its shoulders, anisotropy in the spectral weight, etc., 
can be reproduced in reasonable range of parameter val- 
ues, only by assuming that the system is in the charge 
disproportionated ground state. We have discussed the 
possible reconciliation with experimental data suggesting 
the presence of uniform ladders. 
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Fig. 4. Same as Fig. 2 but for the case with A / 0. Here we 
assume that the sites correspp«<jfejLg to V 5+ ions on the zigzag 
structure of ordered electronst3 E3' have a nonzero value of A. 
Values of the ladder parameters of Set A are used. 



Finally, let us discuss the controversy in the valence 
state of V ions h^Na-V^Os at T > T c : recent^ X-ray 
structural analysescHjo' and NMR measurements© seem 
to suggest the uniform valence of V 4 5+ , the system be- 
ing the coupled uniform ladders with no charge dispro- 
portionation, which apparently contradicts our results 
shown above. However, we immediately notice that, if 
the dynamical fluctuation of the valence state of V ions 
is assumed, a possible way of reconciliation is to take 
into account the difference in the time scale of the mea- 
surements: if the charge (or valence) fluctuation of V 
ions is slow compared with the frequency of electric field 
of light used in the optical measurements, one finds the 
system to be charge-disproportionated, but if the charge 
fluctuation is fast enough compared with the time scales 
of the NMR measurement (an order of 1 — 100 MHz) and 
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